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I. INTRODUCTION 


Ising on fractals. The well known exact solutions of the Ising model in one and two 
dimensions are the only exact solutions we have to date [1, 2]. Ising models on fractals of 
dimension between one and two are natural possible candidates to enter this restricted 
group of solvable models. Actually, we already have solutions on some fractals, those 
with finite ramification number like the Sierpinski gasket, but these are of limited interest 
since they resemble the one dimensional case as they do not possess any phase transition at 
finite temperature [3] . Instead, a fractal with infinite ramification number as the Sierpinski 
carpet, which we know has a non-zero critical temperature [4,5], has to date been studied 
mostly numerically, and few analytical studies are available. In this paper we try to fill 
this gap presenting an analytical study of the Ising model on fractals of dimension below 
two, which include both the gasket and the carpet. We present a method in principle 
able to determine the critical temperatures exactly for all these fractals. Our approach is 
based on approximating the Ising model on non-periodic fractal lattices with a sequence 
of Ising models on periodic lattices. We do this by exploiting our ability to readily solve 
the two dimensional Ising model on an arbitrary periodic lattice using an extension of fhe 
combinatorial method of Feynman-Vodvickenko [6-8]. 

Universality. The understanding of universality classes in dimension equal or above 
two is now quite robust, in particular for system with Z 2 symmetry [9,10]. Less clear is 
the situation in dimension below two and greater than one. Continuous methods [11,12] 
give a fairly good description, and in some dimensions real space renormalization group 
studies are available [13,14], but an explicit solution of the Ising model in some fractal case 
will provide strong indication regarding the reliability or not of continuous methods in 
dimension below two. In fact it is not completely clear if there is a difference between the 
values of critical exponents one can obtain with continuous methods, which usually make 
a continuation of the integer number of dimensions to fractional values, and the actual 
values obtained studying the analogous systems defined directly on lattices of non-integer 
fractal dimension. Another open question is if there is a lower critical dimension for the Z 2 
universality class, or even if this concept is well defined since it might be that universality 
depends on the fine details of the fractal [15]. To clarify all these question, a better 
understanding of the Ising model in fractal dimension, which is a representative of fhe 
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Za universality class, will be of the utmost significance. A further reason why the Ising 
model universality class is interesting in dimension below two is because it is the only 
non-trivial one in the family of the 0{N) models due to the generalised Mermin-Wagner 
theorem [16]. 

Summary of the paper. In Section II we explain how to generalize the combinatorial 
method of Feynman and Vodvickenko to arbitrary periodic lattices. Then, in Section III, 
we explain how to apply it to approximate fractals. After some analytical results, we turn 
to numerical methods to extract the approximate critical temperatures and correlation 
lengths for many different fractals. We finally compare our results with the numerous 
existing Monte Carlo estimates available in the literature and we briefly discuss possible 
future applications of our method. 


11. SOLUTION ON ARBITRARY PERIODIC LATTICES 

A. The model 

Definitions. We briefly review the definitions that specify the model. We consider an 
arbitrary periodic lattice A where at every lattice site i there is a spin variable U; e {-1,1} 
and we define a microstate by a spin configuration {cr}. We assume nearest neighbour 
interactions so that the energy of a given spin configuration is given by 

Ela} = -jY,aiaj. ( 1 ) 

(y) 

If/ > 0 the interaction is ferromagnetic, while it is antiferromagnetic if / < 0. The partition 
function is the sum over all spin configurations weighted by the Boltzmann-Gibbs factor 

= ( 2 ) 

hi {i,j) 

We define f - l/ksT, where ks is the Boltzmann constant, and we set / = 1 since we are 
going to consider only ferromagnetic interactions. 

High temperature expansion. Using the high temperature expansion, the partition func¬ 
tion (2) can be rewritten as 

ZA = 2^(coshff'0A{v), (3) 
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where v = tanh jS and N = Ng is the total number of lattice sites while Ni is the total number 
of links. The function Oa(u) is the generating function of the numbers which count the 
graphs with even vertices of a given length that can be drawn on the lattice A. In this way, 
the problem of solving the Ising model on A is reduced to the combinatorial problem of 
counting even closed graphs on A. In the thermodynamic limit it is the function OA(y) 
which develops the non-analyticity that characterises the continuous phase transition. For 
this reason in the following we will focus on it and disregard the pre-factors appearing 
in (3). 

For high temperature expansion studies of the Ising model on fractals, see [17]. As ex¬ 
plained in the next section, we instead resum the high temperature series by generalizing 
the approach of Feynman-Vodvickenko. 

B. Feynman-Vodvickenko method 

Exact solution on arbitrary periodic lattices. Feynman [7] and Vodvicenko [6] introduced 
a trick to reduce the problem of counting closed graphs to a random walk problem. More 
precisely, the generating function OA(t’) can be computed by counting closed weighted 
random walks paths on A, where the weights are complex amplitudes constructed so that 
the mapping from the high temperature expansion to the random walk problem works 
out correctly [8, 18-20]. 

From the knowledge of the transition matrix Wa of the random walk problem, one can 
determine the explicit form of the generating function from the following relation [8]: 

<I>a{v) = exp J" log det [I - WA(k)]| , (4) 

where the k = {kx,ky) integration is over the region 0 < kx < 2n and 0 < ky < 2n. The 
singular non-trivial part of the free energy for spin logOA for the lattice A, can 

finally be written as 

P/a W = - J J ^ 'os ' <5) 

where we have defined the determinant 

PA(u,k) = det[I-uWA(k)] . (6) 

The matrices WA(k) are m xm matrices with m — sxl, where s is the number of sites in 
the basic tile and / is the total number of links in the basic tile. 
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Critical temperature. If a phase transition takes place, the critical temperature can be 
determined as the real solution of 


Pa{v,0) = 0, 


(7) 


in the range 0 < < 1. Solutions of equation (7) for all Archimedean and Laves lattices 

have been studied in [8]. Equivalently, we can determine the critical v as the inverse of 
the largest positive real eigenvalue of Wa(0). This characterisation is very useful when 
the computation of the characteristic polynomial (7) becomes too demanding. Near the 
critical point, and in terms of the reduced temperature t - T/Tc - 1, the critical exponent 
a is defined by the scaling fj^{t) ~ In particular, a logarithmic singularity of the free 
energy, as the one present in (5), is encoded in a = 0. 

Correlation length. The correlation length can be computed from the knowledge of 
the lattice mass since <^a = I/^a- This last is defined by the following small momenta 
expansion of the determinant. 


Pa{v, k) = Za{v) [m\{v) + k^ + 0{k ^)\, 


( 8 ) 


where Za{v) = ^PA{u,'k)\j^_^ is the wave function renormalization. The lattice mass can 
then be written as m^^{v) = Pa(c’/ Q)IZa{v) and the correlation length takes the form 


Ia{v) = 


j^PA{v,'k)\ 


k=0 


PaM 

The correlation length critical exponent is defined by the relation <^a 
scaling region near the phase transition. 


(9) 


|f| ^ valid in the 


III. FRACTALS 

Approximate solutions on fractals. We now want to use our ability of solving the Ising 
model on an arbitrary periodic lattice to find approximations to the same problem but on 
fractal lattices of fractal dimension below two. 

We will study Sierpinski carpets that are defined in an iterative way. Let us consider 
a two dimensional L x L tile, where some of the squares have been removed: we call 
this the generator of the fractal. It is also the first iteration, A: = 1, of the sequence that 
defines the fractal. Then, given the sequence at iteration k, the next iteration, k + 1, can be 


5 





FIG. 1. Illustration of sc(3, l)k for k-1,2, 3,4. 

constructed by replacing every existing square at the A:* iteration with the generator. In 
the limit k ^ oo this defines a fractal. 

In our approach, the lattice Ak at iteration k is defined as the infinite repetition of the 
tiling obtained at level k. We assume the Ising model on the limiting fractal Aoo = limjt^co Ajt 
exhibits the same behaviour as the Ising model defined on the fractal. In other words, we 
approximate the fractal's determinant PaA'^> 1^) using periodic approximations 


PAjv,k) = limPA,{v,k). (10) 

A:—>oo 

Furthermore we define sc(L, b)i to be the generator where from L X L tile abxb square is 
removed from the center. Then sc(L, b\ denotes the tiling at iteration k. For illustration 
see Figure 1. We will denote with SC(L, b\ the lattice obtained by tessellation of the plane 
with tile sc{L,b)k. We will also study Sierpinski gaskets. Theirs generators are a Lx L 
tile where a single (L - 1) X (L - 1) block has been removed. Sierpinski gaskets have a 
finite ramification number, i.e. one can remove arbitrarily large pieces by cutting a finite 
number of links. 

We denote with Tj^ the critical temperature of the lattice A^ where A is understood from 
the context and similarly for the correlation length 


A. Sierpinski carpets 

Explicit form ofW. We can construct the transition matrix for any finite 

iteration Sierpinski carpet A^ exactly. Explicitly, this can be constructed as the adjacency 
matrix of a weighted directed graph, in which each node of the graph represents one of 
the four directions U,L,D,R associated to each site of the generator of Aj^. In particular, 
when kx = ky = 0 we can represent Wa^(0, 0) as the graph 
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where the indices are modulo L and a = e'? is the complex amplitude required by the 
Feynman-Vodvicenko method. Basically, in terms of directions, clockwise arrows have 
amplitude a, counter clockwise arrows have amplitude a ^ while self connections have 
amplitude one. The weights Wi^j are chosen equal to the matrix representation of the 
generator of the lattice setting Wi^j = 1 if fhe site {i, j) of fhe generator exists and w, ^ = 0 
if it is depleted. The case of the standard Ising model on square lattice can be represented 
as in Figure 2. The momentum dependence of ^y) is obtained by multiplying the 

links outgoing from U with from L with from D with and from R with 

Although we are here interested in approximating fractals, by properly choosing W; y 
this construction gives the transition matrix for any lattice with rectangular tile. For 
example, our construction encompasses the exact solution of the Ising model on all 
possible two dimensional depleted lattices, including the case of a random basic tile. 

Standard Ising model. We shortly review the solution of the standard two dimensional 
Ising model to exemplify the method. From Figure 2 we immediately reconstruct the 
transition matrix 


^^Ising{kxr l^y) — 



Ipiky 

0 

ae4 


gtkx 

Ipikx 

0 

0 

ae~4 

e-iky 

X. 


0 


g-ikx 


( 11 ) 


with a = . The determinant is readily computed and gives the well known Onsager's 

solution [2]: 


Plsingiv, k) = (l + - 2 U (l - (cos kx + cos kyj . 


( 12 ) 


Setting Pisi„g{v,0) - {1 — 2v - tPY = 0 gives Vc — 0.414214... as the only solution in the 
range 0 < u < 1; this correspond to the critical temperature = 2.26919... first computed 
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FIG. 2. Graph representing the transition matrix WisingiO, 0) for the standard Ising model on a 
square lattice. 


by Kramer and Wannier using duality arguments [21]. Finally, using (9) we find the exact 
form for the correlation length 

^lsing{v)^l - , , (13) 

which clearly diverges for T ~ since the denominator vanishes. 

To check our formalism we can try to solve a redundant version of the two dimensional 
Ising model, defined on a tile of size L > 1. With our previously defined notations, they are 
equivalent to SC(L, 0)i. For example in the cases L = 2,3 we find to following characteristic 
polynomials, 

^’sc( 2 , 0 )i(t^, 0 ) = (1 + v^f{l + lv- i;^)^(l -Iv- v^f 

Tsc( 3 . 0 )i {v, 0) = (1 - 2z; - v^f{l + lv + 2v^ - 2v^ + v^il -v + 2v^ + v^ + , (14) 

which indeed have only the = 0.414214... solution in the range 0 < z; < 1. This represents 
a non-trivial check of our ability to construct the transition matrices for an arbitrary lattice. 
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The momentum dependence becomes rapidly very complicated but a similar analysis can 
be made for the correlation length. 


B. Analytical results 


Analytical solutions for the Sierpinski gaskets. Analytical relations can be found for all 
the Sierpinski gaskets defined on a L x L grid. The cases L = 2,3,4,5,... are shown in the 
legend of Figure 3. We are able to give the analytical form for the determinant: 

Pa^'u, k) = - 2 (cos kx + cos ky^. (15) 


Sincey < 1 the infinite iteration limit leads just to one lim)t^c>oFA^(z’,k) = 1. ThusOA^(y) = 1 
as in the one dimensional case and the singular non-trivial part of the free energy per spin 
is zero. We recover in this way the result that Ising models on Sierpinski gaskets do not 
magnetise [3, 22]. 

Even if Too = 0, it is interesting to infer from (15) the exact critical temperature for any 
finite k: 


U(V2-1)‘'‘* ' 


( 16 ) 


This relation is instructive since it shows that convergence to the limiting value is very 
slow, more precisely logarithmic, as can be seen in Figure 3. It is interesting to look also 
at the correlation length, which from Eq. (9) turns out to be 


Uv) = 2 






(17) 


This relation is visualised in Figure 4. The correlation length critical exponent is one, as 
in the one dimensional case, for all Sierpinski gaskets. 

It is clear that a similar analysis, with similar conclusions, can be made for other 
families of fractals with finite ramification number. It is probably possible to obtain a 
closed formula for Pkiv, k) for any fractal with this property. This is a clear indication of 
their triviality and effective one dimensional behaviour. 
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FIG. 3. The critical temperature for the Sierpinski gaskets as a function of k. Convergence towards 
Too = 0 is very slow (logarithmic). The k = 3 generators for (from left) the L = 2,3,4,5 Sierpinski 
gaskets considered in the study are shown in the legend. 


Analytical solution of thek = 1 Sierpinski carpet. We can give the analytical solution of 
the Sierpinski carpet SC(3, l)i. The critical temperature is the solution of 


Psc{ 3 ,ih 0) = (l - - 16v^ - lOv^ - 20v^ + u^- 24i;^ + . (18) 


The only root in the range 0 < u < 1 is = 0.495968... which gives = 1.83842... as 
reported in the k = 1 entry of Table I. Note that it is a non-trivial fact and a consistency 
check that the 12^^ degree polynomial in Eq. (18) has only one real solution in the physical 
range. 

In this case we are also able to determine the full momentum dependence of the 
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FIG. 4. The correlation length for the L = 3 Sierpinski gasket as a function of T for the values 
= 1,2,10,50. 

determinant 

-PsC(3,l)i (^/k) = 

1 + lOv^ + 20i;^ + 119i;^ + 324v^^ + 876i;^^ + 1284z;^^ + 983u^^ + 412i;^^ + 58v^^ + 8i;^^ + 
-4u^(l - vY{l + 6v^ + 21v^ + 52v^ + 69v^ + 72v'^° + 29v^^ + 6t;^^)(cos 3k^ + cos 3ky) 
-2v^{l - (7 + 18v^ + 24z;^ + 14z;^ + (cos 3{kx + ky) + cos 3{kx - ky)^ 

+2v^(l - (1 + 4v^ + 3v‘^^ (cos 6kx + cos 6ky^. (19) 

This relation clearly illustrates how non-trivial are the explicit solutions already at the level 
of the first iteration. It also shows how higher harmonics are excited, and that in Pa^{v, k) 
the coefficients of the trigonometric functions are polynomials in v. We have obtained 
similar relations for many k — 1 non-trivial fractals, i.e. with infinite ramification number, 
while we have not been able to obtain closed analytical forms for the determinants as a 
function of k, and we suspect this to be a formidable task, even if not hopeless. Such a 
closed formula will constitute an explicit exact solution of the model. 
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FIG. 5. The L = 3 Sierpinski carpet(s) considered in the study at their fc = 4 approximation. 


Finally, we also report the correlation length in the k — 1 case. Inserting Eq. (19) into 
Eq. (9) gives 


4lW = 2 


+ 12!)“ + 18d“ + 58d“ - 13!)“ + 144!)“ - 20!)® + 138!)® + 7!)^ + 104!)® + 2!)® + 42!)'* + 5 jfi + \2v^+ 2) 
f !)12 + 2 !)!® - 24 !)^ + !)® - 20!!2 - 10 !)® - 16 !)® + 5 !)^ - 4 !)® + 


1 

2 


This correlation length diverges consistently at = 0.495968... and when expressed in 
terms of the reduced temperature is plotted as the upper curve in Eigure 9. Clearly v = 1 
as expected from universality. 


C. Critical temperatures 

Numerical analysis of critical temperatures. We have reduced the calculation of the crit¬ 
ical temperature on a lattice Ak to finding the largest positive real eigenvalue of the 
matrix Wa^( 0 , 0 ) corresponding to the weighted adjacency graph defined in the previous 
section. The X 4L^*^ matrix is sparse, but its size grows rapidly as a function of 
the k. Eor k > 1 we are not able to calculate analytically the eigenvalues but we instead 
resort to numerical calculations. 

We use the shifted Arnoldi solver of Mathematica, which uses the ARPACK library 
and is sufficient for an initial proof of concept. The algorithm can be used to compute 
an arbitrary number of eigenvalues in the neighbourhood of a complex number, usually 
referred as the shift parameter. In our study, we compute the eigenvalue Aa^ using as a 
shift the eigenvalue Aaj^.j • 

A caveat of our approach is that it would give a wrong estimate of Aa^ if there was a 
complex eigenvalue with a small imaginary part in the neighborhood of Aa^_i • However, 
it turns out that the eigenvalue is isolated, as for instance illustrated in Eigure 6 in the case 
of SC(3, 1 ) 3 , for which we can compute the entire spectrum of the matrix Wsc( 3 ,i) 3 ( 0 , 0 ). 
Eor larger k we investigate the stability of our prediction depending on the value of the 
chosen shift parameter. Eor all the cases considered we observed that changing the shift 
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FIG. 6. Eigenvalues of Wsc( 3 ,i )3 • The only real eigenvalue larger than one is shown using a blue 
diamond. 

leads to the same estimate of the critical temperature when we are able to compute enough 
eigenvalues. If we could not compute enough eigenvalues, then we find only non real 
eigenvalues. Therefore, in practice this numerical limitation does not arise. 


Using this procedure, we can calculate the critical temperature up to A: = 7, for the L = 3 
fractals. To achieve this, we need to find a specific eigenvalue of a 19 131 876 x 19 131 876 
matrix. The calculations are limited by the available memory. All the computations 
performed in this section have been achieved running Mathematica on a single node with 
20 cores and 128GB of memory. Without further method improvements a machine with 
more memory would be needed to compute Tk for k>7. 

The results for the L — 3 Sierpinski carpets are given in Table I. They are illustrated 
for the 8 fractals considered in Figure 7. The fractals considered here have two different 
fractal dimensions. However they differ by there number of active bonds. The three 
fractals that have a dimension close to two show a much faster convergence than the 
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FIG. 7. Exact critical temperatures for the first seven approximands to the non-trivial L = 3 
Sierpinski carpets considered, together with a fit. 

others. As seen in the previous section, in the case of the Sierpinski gaskets, fractals can 
show a very slow convergence in the A: —> oo limit. We thus see that, for the fractals 
considered, the smaller the fractal dimension, the slower is the pace of convergence. We 
do not attempt to perform any extrapolation to A: —> oo since, in our set-up, we lack a 
scaling theory that governs this limit as the volume is already infinite for every k. 

As mentioned earlier, the method can be easily applied to other fractals, and we 
investigate some of them with L — 4,5 and 7 as summarized in Table II. The fractal 
SC(7,3) is defined by removing nine distinct uniformly distributed cells as illustrated in 
Figure 8. The two fractals with L = 7 generators considered here have the same fractal 
dimension but different lacunarity. They have been discuss in [15, 23]. 
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Generator 

df 

A = 1 

(N 

II 

CO 

II 

II 

A = 5 

kO 

II 

13^ 

II 

m 

1.893 

1.83842 

1.67971 

1.61601 

1.58935 

1.57798 

1.57310 

1.57099 

a 

1.893 

1.83842 

1.66680 

1.59188 

1.55769 

1.54140 

1.53319 

1.52872 

□ 

1.893 

1.83842 

1.65386 

1.56759 

1.52566 

1.50446 

1.49331 

1.48719 

X 

1.771 

1.48866 

1.18962 

1.04440 

0.965875 

0.920115 

0.892608 

0.875999 

a 

1.771 

1.48866 

1.18310 

1.03567 

0.955384 

0.908286 

0.879960 

0.862996 

# 

1.771 

1.29944 

0.983021 

0.830078 

0.739657 

0.679436 

0.636087 

0.603146 


1.771 

1.29944 

0.97052 

0.80394 

0.699862 

0.626109 

0.569437 

0.523428 

d 

1.771 

1.29944 

0.958433 

0.780739 

0.667582 

0.586519 

0.524012 

0.473526 


TABLE I. Exact critical temperatures for the various iterations of the L - 3 Sierpinski carpets. The 
fractals are ordered from the largest to the lowest critical temperature. Note that for k = 1 different 
lattices have the same critical temperature, since at this level, there are exactly the same. 

Fractal df k = l k = 2 k = 3 fc = 4 k = 5 
SC(4,2) 1.792 1.62129 1.36891 1.25015 1.18451 1.14280 
SC(5,1) 1.975 2.11926 2.07899 2.06904 2.06672 
SC(5,3) 1.723 1.48748 1.19857 1.05787 0.97483 
SC(7,3) 1.896 1.92863 1.85117 1.8334 1.82927 
SC(7,3) 1.896 1.57100 1.39728 1.34601 1.32719 

TABLE 11. Exact critical temperatures for the various Sierpinski carpets with L > 3. 

Comparison with Monte Carlo approach. Defining Ajt,/ = sc(L, b)k,i as a finite lattice de¬ 
fined by an / X / array of sc(L, h\ building blocks. Defining Tk,i as the critical temperature 
of the lattice Ajt,// where obviously for a finite I the critical temperature is defined for 
instance as the maximum of the specific heat. Usually Monte Carlo studies reports Tjt,! 
whereas we compute Tjt.oo = lim;^oo While Tjt,; depends on the lattice definition of the 
critical temperature, the value of oo is unique. The critical temperature of the fractal is 
approached in the limit of A: ^ oo in both cases. Since it is proven that for fractals with 
infinite ramification number Too > 0, the two approaches must yield to the same limiting 
value. The rate of convergence is a priori unknown in both cases. 

In table III we compare our results for Tk,oo with the results for Tk,i obtained using Monte 
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FIG. 8. Generators sc(7,3)2 (left) and sc(7,3)2 (right). 

Carlo simulations for various fractals. The lattice estimates of the critical temperature are 
in a good agreement with our results and almost always within the estimated errors. As 
expected from the previous considerations, the fractals with the highest fractal dimensions 
exhibit also a better agreement. 


D. Correlation lengths 

As explained in section IIB, we can also estimate the correlation length of fhe system 
as a function of k,t — T/Tq-I (the reduced temperature) and Ajt using Eq. (9). We expect 
from universality that all approximands have v = 1, since for finite k, they all belong to 
the universality class of the two dimensional Ising model. Different critical exponents for 
the fractal can emerge only in the limit k ^ oo where the type of singularity manifested 
by the free energy at the critical point can change. 

We checked numerically this expectation. Evaluating Eq. (9) exactly for various values 
of the reduced temperature, we computed the correlation length for various L - 3 fractals 
with k up to four. Our results for SC(3,1)^ are illustrated in Eigure 9 where the correlation 
lengths have been normalized to exhibit universality. As can be seen all the curves are 
compatible with v = 1, but the scaling region shrinks as k increases. Our results for fhe 
normalized correlation lengths of all L = 3 fractals considered in Table III are represented 
in Eigure 10. This picture represents a strong confirmation of universality. 

This analysis implies that if the limit k ^ oo is continuous then the critical exponent v 
is one for the fractals of dimension below two. In contrary, if the limit is discontinuous 
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Authors 

Tc 

V 

k 

SC(3,1) df 

= 1.8927 



Bonnier et al. (1987) [26] 

1.54 

1.3 

3 

Our work 

1.56759 

- 

3 

Pruessner et al. (2001) [25] 

1.5266(11) 

- 

4 

Our work 

1.525660 

- 

4 

Pruessner et al. (2001) [25] 

1.5081(12) 

- 

5 

Our work 

1.504460 

- 

5 

Pruessner et al. (2001) [25] 

1.4992(11) 

- 

6 

Bab et al. (2005) [27] 

1.4945(50) 

~ 1.39 

6 

Our work 

1.493310 

- 

6 

Carmona et al. (1998) [28] 

1.481 

1.70(1) 

7 

Monceau et al. (1998) [29] 

1.482(15) 

1.565(10) 

7 

Our work 

1.48719 

- 

7 

Monceau et al. (2001) [30] 

1.4795(5) 

> 1.565 

8 


SC(4,2) 

df 

= 1.7925 



Carmona et al. (1998) [28] 

1.077 

3.23(8) 

6 

Monceau et al. (2001) [30] 

< 1.049 

>3.37 

6 

Our work 


1.1428 

- 

5 

SC(5,1) 

df 

= 1.9746 



Monceau et al. (2001) [30] 

2.0660(15) 1.083(3) 

5 

Our work 


2.06672 

- 

4 

SC (5,3) 

df 

= 1.7227 



Monceau et al. (2001) [30] 

< 0.808 

>4.06 

5 

Our work 


0.974828 

- 

4 

Ising 2d 


2.269 

1 

- 


TABLE III. Comparison with the literature (Monte Carlo study). 
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FIG. 9. Normalized correlation lengths of the Sierpinski carpet starting from above for k = 1,2,3,4 
as a function of the reduced temperature t. As expected from universality all all curves tend to 
for small enough t. 

other values are possible, but our approach cannot determine them, unless we are able 
to calculate analytically the determinant PA^iv,'k.) as we did for the Sierpinski gaskets. 
Without such an analytical formula, at any finite k, it is difficult to estimate the critical 
exponent v of the fractals. 

Comparison with the Monte Carlo approach. Defining, in analogy with what done in the 
case of the critical temperatures, Vk,i as the critical exponent for the iterations of systems 
of size 1. We compute Vk,co, and find Vk,oo = 1 for all k as expected from universality. For 
the fractal critical exponent Voo = limjt^co Vk,i to be different from one the limit has to be 
discontinuous. The standard theory of finite size scaling [24] applies to changing / while 
keeping k fixed, and this scaling should lead to lim/^oo W; = 1 as required by universality. 
This was already noted by [25]. To our knowledge there is no scaling theory with respect 
to k at fixed I, in particular for / = 1, which is used in some Monte Carlo simulations. For 
/ = oo there is not such theory because Vk,oo = 1 for all finite k. Hence, some Monte Carlo 
simulations rely on a possible scaling on Vjt,!/ which might not exist. 

Furthermore, our results show that the scaling region where ^ ~ t~^ shrinks as k is 
increased, as shown in the Figure 9. This suggests that it becomes increasingly difficult 
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FIG. 10. Universality strikes back: normalized correlation lengths of the non-trivialL = SSierpinski 
carpets for fc = 4 as a function of the reduced temperature t. As expected from universality all 
curves tend to (dashed lines) for small enough f. 

to compute the critical exponent i keeping k fixed and calculating it using the reduced 
temperature as the scaling variable, which is an other approach used in Monte Carlo 
simulations. 

As shown in the table III, Monte Carlo simulations provide estimates for the critical 
exponents, in particular v. They report values of v > 1 even bigger than four [30]. It 
might be that the lattice simulations are able to capture the right universal properties of 
the fractals. Nevertheless, a better theoretical understanding of the situation is needed. 


IV. CONCLUSIONS AND OUTLOOK 

We showed that it is possible to approximate the solution of fractal Ising models via a 
sequence of exact solutions of Ising models on finite periodic representations of the fractal 
under consideration. We found that the rate of convergence to the exact solution can be 
very slow, as the explicit example of the exact solution of the Sierpinski gasket model 
shows. We found estimates for all L = 3 fractals, and in particular for the Sierpinsky 
carpet. Numerical improvements can rapidly refine our results ultimately leading to 
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the accurate determination of the exact critical temperatures for all non-trivial fractals of 
dimension smaller than two. 

The problem is much more difficult in the case of the critical exponents since universal¬ 
ity ultimately sets in and renders any finite periodic approximation useless to the scope. 
But we can still speculate on the actual values for the critical exponents of the fractals. 
Our analysis suggests two scenarios, related to the fate of the limits Vco = limjt^oo Vk,oc and 
aco = limjt^oo «fcoo. Assuming the hyper scaling relation a = 2-dfV to hold true for fractals 
with infinite ramification number, then the continuity of the first limit implies Voo = 1 and 
a change in the singularity structure of the free energy, i.e. ctoo -2 - df. In this case both 
critical exponents will be continuous in the df ^ 1 limit. If, instead, there is no change of 
singularity structure, i.e. Poo(u,k) retains the polynomial zero it has in the finite k cases, 
then aoo = 0 and, again assuming the hyperscaling relation, we find Voo = 2/d/, which 
is not continuous in the d/ —> 1 limit but is a fairly good approximation for all fractional 
dimensions as compared to the RG results [11, 12]. Obviously, it can also be that both 
limits are discontinuous or that the hyperscaling relation is either violated or does not 
contain df. It can also be that these limits are sensible to the type ot fractal under study, i.e 
they depend on fine details such as lacunarity or cormectivity [15, 31]. 

The question of the existence of a lower critical dimension, for fractals with infinite 
ramification number (or more restrictive properties), can now in principle be addressed 
by our method once the numerical routines are improved, since we have seen that the rate 
of convergence of the critical temperatures Tk is slower the smaller the fractal dimension 
is. High values of k will thus be needed to resolve the neighbourhood of a possible lower 
critical dimension. 

Further interesting applications of our method are related to the study of Ising models 
on other non-translationally invariant lattices, like those defined on aperiodic or random 
lattices, or with random interactions. 
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